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ABSTRACT 

We analyze 1^-band photometry of the aperiodic variability in T CrB. By applying a 
simple idea of angular momentum transport in the accretion disc, we have developed a 
method to simulate the statistical distribution of flare durations with the assumption 
that the aperiodic variability is produced by turbulent elements in the disc. Both 
cumulative histograms with Kolmogorov-Smirnov tests, and power density spectra 
are used to compare the observed data and simulations. The input parameters of the 
model R[ n and a are correlated on a certain interval and the most probable values are 
an inner disc radius of R m ~ 4 x 10 9 cm and a viscosity of a ~ 0.9. The disc is then 
weakly truncated. We find that the majority of turbulent events producing flickering 
activity are concentrated in the inner parts of the accretion disc. 
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1 INTRODUCTION 

Cataclysmic variables (CVs) are interacting binaries with a 
white dwarf as the primary component and a late type main 
sequence star as the secondary star which fills its Roche- 
lobe (see Warner 1995 for review). Consequently mass is 
transferred to the white dwarf through the inner Lagrangian 
point (Li) and an accretion disc is formed (not in the case of 
magnetic polars) . The accretion disc is the site where several 
known physical mechanisms manifest themselves in the form 
of short-term stochastic variations in brightness, generally 
known as flickering. These oscillations are one of the most 
characteristic observational properties of accreting systems 
with typical amplitudes ranging from a few dozen of milli- 
magnitudes up to more than one magnitude on timescales 
from seconds to minutes. Symbiotic stars are related systems 
to CVs where the secondary star is an evolved red giant 
where the mass loss can be due to either the wind from the 
giant star or the Roche-lobe secondary star, as in "classical" 
CVs. Thus accretion discs can be formed with the associated 
stochastic variations. 

The first systematic study of mechanisms generating the 
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flickering was performed by Bruch (1992), where he com- 
pared theoretically estimated energies and timescales for in- 
dividual mechanisms, with observationally derived values for 
a selection of CVs and related symbiotic stars. The author 
proposed four possible mechanisms responsible for the ob- 
served variations in brightness: (1) unstable mass transfer 
from the Li point and interaction with the disc edge, (2) dis- 
sipation of magnetic loops, (3) turbulences in the accretion 
disc and (4) unstable mass accretion onto the white dwarf. 
Apparently it was not enough to restrict the interpretation 
of the flickering in the whole family of CVs to only these 
four mechanisms. For example, in the case of a truncated 
accretion disc (intermediate polars) or in the case of the to- 
tal absence of a disc (polars) the source of such variations is 
known to be the accretion stream threaded by the magnetic 
field and its subsequent interaction with the surface of the 
white dwarf - accretion shock. 

As mentioned earlier, Bruch (1992) performed a sys- 
tematic study of flickering activity in a sample of CVs and 
symbiotic systems. The promising models which do not con- 
tradict the empirical results, are unstable mass accretion 
onto the central accretor or turbulence in the inner accre- 
tion disc. Bruch (1996) studied the variations of photomet- 
ric data as a function of orbital phase in Z Cam. The source 
of the flickering was localized in the vicinity of the central 
accretor, but the outer disc edge and its interaction with 
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Table 1. The observing log. T is the duration of the observation 
in hours and N is the number of data points in the lightcurve. 
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Date 


HJD-2450000 


T 


N 






(start) 


[h] 




1 


08.03.1996 


150.51153 


1.94 


120 


2 


09.03.1996 


151.44714 


1.80 


119 


3 


19.04.1996 


192.35582 


3.85 


240 


4 


20.04.1996 


193.38731 


3.74 


240 


5 


22.04.1996 


195.48539 


2.30 


150 


6 


23.04.1996 


196.34682 


2.79 


179 


7 


03.07.1996 


267.37947 


1.29 


89 


8 


11.12.1996 


428.61449 


1.80 


120 


9 


15.01.1997 


463.54001 


3.87 


219 



the gas stream from the Li point cannot be neglected (hot- 
spot). Bruch (2000) studied 4 eclipsing systems: HTCas, 
V2051 Oph, UXUMa and IP Peg. In these systems, similar 
to Z Cam, the flickering source is located in the central part 
of the disc or the hot-spot. For V2051 Oph the region of the 
flickering was confirmed by Baptista and Bortoletto (2004) 
using eclipse mapping. The hot-spot is the source of low fre- 
quency flickering and the central disc of the high frequency 
flickering. Finally, T CrB case was studied by Zamanov and 
Bruch (1998). The U photometric data showed fast vari- 
ability located near the central white dwarf. Following these 
observational results it is clear that the dominant source of 
flickering variability in CVs and TCrB is the central part 
of the disc or the vicinity of the white dwarf. A promis- 
ing scenario to explain this phenomenon is the existence of 
magneto-hydrodynamic turbulence transporting the angular 
momentum outwards (Balbus and Hawley 1998). However 
why only the central part of the disc is the source of the 
observed turbulence? 

In this paper we focus on the analysis of turbulent trans- 
port of angular momentum in the disc. We have decided to 
study T CrB due to its rich phenomenology of rapid variabil- 
ity (i.e. flickering activity) and because of it's disc properties 
are suitable for our model. In Sec. [2] we present our observa- 
tions. Our treatment of turbulence in the disc as the source 
of flickering is explained in Section [3] and the simulations 
are described in Section [4] Finally the results are analyzed 
in Sec. [5] and discussed in Sec. [6] 



2 OBSERVATIONS 

The observations of TCrB were taken with the 60 cm 
Cassegrain telescope at the Skalnate Pleso observatory of the 
Astronomical institute of Slovak Academy of Sciences. We 
used a single channel photoelectric photometer with UBVR 
Johnson filters. We continuously cycled the UBVR filters 
with 10 s integrations per filter and 1 s dead time for each 
filter change, which results in a time resolution of 45 s per 
filter. The data were first used in the long-term lightcurve 
analysis by Hric et al. (1998). 

Figure Q] shows two example runs taken in the Johnson 
V filter. The individual lightcurves consist of rapid aperiodic 
variability with durations of ~ 10 2 — 10 3 seconds. We have in 
total ~ 24 hours of fast photometry divided into 9 observing 
runs. Tab. [T] summarizes all our observations. 
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Figure 1. Two examples of our V-band lightcurves. The upper 
panel shows data taken on the 23 April 1998 and the lower panel 
shows datas taken on 15 January 1997 

3 SOURCE OF APERIODIC VARIABILITY 

T CrB is a non magnetic system, therefore the source of the 
aperiodic variability must be attributed to one of the four 
mechanisms proposed by Bruch (1992). Using integration 
times of 10-15 minutes, Zamanov et al. (2005) discovered 
variability in the H a emission line which has a character- 
istic timescale of ~ 1 hour. Using the velocity information 
the authors localize the variability at ~ 20 — 30 Rq. This 
is consistent with the disc radius and the gas stream from 
the L\ point can be the source of this variations. Our data 
also exhibits variability but on timescales shorter than ~ 25 
minutes and therefore must have a different origin. Further- 
more Zamanov and Bruch (1998) and Zamanov et al. (2004) 
claim that flickering in TCrB with the same timescale as 
ours arises from the vicinity of the white dwarf. Therefore, 
we conclude that unstable mass transfer from the Li point 
is not likely to be the source of our observed fast aperi- 
odic variability. In our paper we focus on the transport of 
angular momentum in the disc. This can be produced by 
magneto rotational instabilities (Balbus and Hawley 1998) 
and we propose that the source of the aperiodic variability 
is associated with magnetohydrodynamic turbulence. 

4 SIMULATIONS 

4.1 The idea of a turbulent element 

Our idea is based on a simple approximation of mass transfer 
in the disc. A turbulent eddy in the disc can be represented 
by a blob of matter with mass m, Keplerian tangential ve- 
locity at distance r 

/GMA 1 ' 2 

Vt ( r )= r {—) (!) 

and dimension scale x. Mi is the white dwarf mass and G 
is the Newtonian gravitational constant. The largest possi- 
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ble dimension scale x of such a turbulent element at radius r 
from the central object is the scale height of the disc H (r) at 
this radius. The blob (which has dimension x and a spher- 
ical shape approximation) penetrates inwards or outwards 
(a blob with excess of angular momentum moves outwards, 
whereas a blob with deficit moves inwards) with a distance 
proportional to x and a radial viscous velocity v T (r). The 
mass of the eddy mt is therefore the product of the volume 
and the density p(r) ~ E(r)/i/(r), where £(r) is the surface 
density at distance r i.e. £(r) = J p(r)dz. Hence, 



4 3 
mt{r) = -ax p(r). 



(2) 



Changing the blob's position from r to r + i or r - x also 
changes it's Keplerian velocity v t and angular momentum 
which is defined by 



L = r x mv t . 



(3) 



Thus by summing up all the turbulent elements we obtain 
the global transport of angular momentum in the accretion 
disc. Every radial penetration of a turbulent element pro- 
duces energy (by "friction" ) which is liberated in the form of 
an observed flare. Therefore, it is expected that the timescale 
t of the penetration is equal to the observed flare timescale. 



4.2 Statistical distribution of the eddy dimension 
scale 

It is assumed that the number of eddies decreases with their 
dimension scale x. Therefore, we can approximate the dis- 
tribution function of eddy sizes f(x) by an exponential func- 
tion (in analogy to the stochastic movement of particles in 
Brownian motion) 



f(x) — A expBx. 



(4) 



where B < 0. The limits of the distribution are /(0) = 1 and 
f(H) — 0.01, where H is the scale height of the disc. The 
minimum of the distribution is not because the exponential 
function Q cannot reach and thus we obtain 



f(x) = exp 



In f(H) 
H 



(•») 



Consequently, the distribution of flare timescales t will be 
given by f(t) with t ~ x/v T (r). 

4.3 Application of the distribution function 

We can define two adjacent rings in a disc, with radius r 
and r + Ar where Ar is the ring radial thickness. The ring 
then has a mass 



m(r) = 27rrArE(r) 
and 

m(r + Ar) = 2-k(t + Ar)ArE(r + Ar) 



(6) 



(7) 



The difference in angular momentum between these two 
rings is AL(r) = L(r + Ar) — L(r), where the angular mo- 
mentum of the adjacent rings is 



L(r) = 2-Rr i T l (r)Ar 
and 



( GM\ 



1/2 



(8) 



L(r + Ar) = 2-k[t + Ar) £(r + Ar)Ar 



GMi 
(r + Ar) 3 



1/2 



(9) 



This difference is produced by discrete blobs of diameter x 
and distribution f(x), moving stochastically in radial direc- 
tions between the two adjacent rings. The angular momen- 
tum of an individual blob of size x located at radius r is 
then given by eqs. Q]-[3] i.e. 

T I \ 4 3 , v 2 / GM\ \ 1 / 2 

L t (r,x) = -iTX p(r)r I —3- I (10) 

whereas for a blob at distance r+Ar, the angular momentum 
is given by 

1/2 



L t (r+Ar,x) — —irx p(r+Ar) (r+Ar) 



GMi 



(r + Ar) 3 



(11) 



Finally the difference in angular momentum is ALt(r, x) = 
L t (r + Ar, x) - L t {r,x). 

Until now the distribution function has a maximum 
/(0) = 1. This would be valid for the whole disc only if 
AL(r) and H(r) were the same (for all r). This is of course 
not true because of the different disc scale height H (r) of 
each ring which determines the largest dimension scale of the 
turbulent element. Consequently f(x) is then also a function 
of r, i.e. f(r,x). The angular momentum difference is differ- 
ent for every pair of rings, not only because of H(r), but 
also because of v t (r). Therefore, f(r,x) must be multiplied 
by some scale parameter k — k(r) which depends on radial 
position of the pair of rings. This parameter is needed for the 
final construction of the flare duration histogram; for every t 
we use k(r)f(r,x/v r (r)) between every pair of rings starting 
with the inner disc radius and ending with the outer disc ra- 
dius (the step Ar must be chosen). Of course, the condition 
t ^ H(r)/v T (r) must be fulfilled, i.e. the flare duration can- 
not be longer than that allowed by the biggest turbulence 
dimension scale. 

The scaling parameter k(r) can be obtained through 
the following simple argument. We can compare AL(r) with 
the sum of AL t (r, x) between one pair of rings at distance r. 
Summing up angular momentum of all the eddies is equiv- 
alent to integrating AL t (r,x) with the distribution f(r,x) 
(equation [5} with the maximum value at 1. We thus obtain 



H(r) 



Al(r) 



f(r, x)ALt(r, x)dx. 



(12) 



where this summed angular momentum difference should be 
equal to the global angular momentum difference AL(r). 
The scaling factor fc(r) is then given by 

k(r)Al{r) = AL(r) (13) 

which gives 

AL(r) 



k(r) = 



Al(r) 



(14) 



Finally, we obtain the distribution of relative duration by 
summing up all durations t for each ring, i.e. 



C(t) = > Mn)f(n,t). 



(15) 



where n is the number of imaginary rings (derived from R ln , 
Rout and Ar) and n are the ring diameters. 
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4.4 The model and free parameters 

For the accretion disc we use the standard Shakura and Sun- 
yaev (1973) model (see e.g. Frank et al. 1992). The disc scale 
height, surface density, volume density and radial viscous 
velocity are then given by the equations 

H = 1.7 W s a-^M^%M-rr% 8 f^ [cm] (16) 
E = ^^^/r [g cm" 2 ] (17) 
p = 3.1 10- 8 a- 7 / 10 Mi 1 c / j°6< 8 '-ro 15/8 / 11/5 [g cm" 3 ] (18) 

„„ ln 4 4/5 ji>3/10 ,,-1/4 -1/4,-14/5 r -ll /in\ 

v r = 2.7 10 a M^-LgMj/ r 10 / [cms ] (19) 
where 



in terms of no = r/(10 10 cm), M 1>s = Mi/(1M Q ) and 
Afacc,16 = Macc/(10 16 g S" 1 ). In our modelling the free pa- 
rameters are the mass of the central star Mi, the viscosity 
parameter a, the inner disc radius Ri n , the outer disc radius 
Rout, and the mass transfer rate from the secondary star 
which equals the mass accretion rate M acc through the disc 
for a steady state disc. The radius of the white dwarf Ri is 
approximated by the formula given in Nauenberg (1972). 



5 ANALYSIS 

5.1 Input parameters 

The inner disc radius is a free (searched) parameter and 
for the starting value we use the radius of the white dwarf 
calculated from Nauenberg (1972). The outer disc radius 
is first approximated by ~ 0.5 of the primary Roche lobe 
radius 7?li calculated from Paczynski (1971) 

* L1 = a462ffl {mTTmT ^ < 21 > 

where a is the binary orbital separation and Mi and Mi are 
the primary and secondary masses respectively. With the or- 
bital period of 227.5 days (Kraft 1958, Paczynski 1965) this 
gives a disc radius of up to ~ 10 12 cm. Using the Shakura and 
Shunyaev (1973) model, for the disc annuli with diameter 
10 12 cm we obtain an effective temperature of ~ 1000 K and 
for a disc diameter of 10 11 cm of ~ 5000 K (see e.g. Frank et 
al. 1992). For effective temperatures less than 5000 K (an- 
nuli with radius larger than 10 11 cm) the contribution to 
the optical emission is not significant and can be neglected. 
Therefore, we can adopt as _R out = 10 11 cm in our analysis for 
the V filter data. 

The white dwarf mass in TCrB is 1.37 ± 0.13 M 
(Stanishev et al. 2004). Using UV observations the authors 
also derive a mass accretion rate onto the central accre- 
tor of Mi n = 2.5 x 10 17 g s" 1 during the low-state and 
M in = 1.9 x 10 18 g s" 1 during the high-state. Accord- 
ing to Fig. 4 in Stanishev et al. (2004), all our observa- 
tions were taken during high-state and therefore, we assume 
Mi n = 1.9 x 10 18 g s _1 . According to the disk instability 
model (see Lasota 2001 for review) this value is higher than 
the critical mass accretion rate for the cold disc (Fig. [2} 
for disc radius less than 10 11 cm, but lower than the critical 



mass transfer rate for a hot disc for larger radius. Stanishev 
et al. (2004) also proposed that the disc is hot and optically 
thick up to a disc radius of ~ 1 Rq ~ 7 x 10 10 cm. But if 
there is an intersection of the mass accretion rate with the 
critical values, the disc is unstable and should exhibit dwarf 
novae activity and cannot be in a steady state. This is not 
observed and all the luminosity variations on timescale of 
days to weeks can be explained by the activity of the red 
giant atmosphere (Zamanov and Bruch 1998). 

Irradiation by the central hot source can stabilize the 
disc on the hot branch (see Lasota 2001 for review) to larger 
disc radii as allowed by the critical mass accretion rate. The 
white dwarf in the T CrB system can have a temperature of 
about 10 5 K (see e.g. Hack et al. 1993) which is very high 
compared to the typical white dwarf temperature in classical 
dwarf novae. Hameury et al. (1999) studied the influence 
of the white dwarf irradiation into dwarf novae cycle. The 
irradiation temperature is given by 

TL = (1 - f3)T?- [arcsinp - p(l - p 2 ) 1/2 ] (22) 

where p = Ri/r, Ti is the white dwarf temperature and 
(1 — j3) is the fraction of the incident flux which is absorbed 
in optically thick regions, thermalized and reemited as pho- 
tospheric radiation (see e.g. Friedjung 1985, Smak 1989). 
Taking the absorbed radiation to be of about ~ 10 % and 
a white dwarf radius of 0.2 x 10 9 cm we get temperature 
of 7.700 K for a distance of 10 12 cm. This temperature lies 
in the region of hydrogen ionisation. Therefore permanent 
hot state is possible up to a radius of 10 12 cm in the case 
of ~ 90 % reflection of incident radiation from the central 
massive and hot white dwarf. 

Finally, Zamanov and Bruch (1998) has claimed that 
M2 = Afin in T CrB, and the disc can be assumed to be 
in a hot steady state. We thus assume that the material is 
hot, ionized and optically thick and the conditions for the 
standard Shakura and Shunyaev (1973) model are satisfied. 
Therefore, we adopt M acc = 1.9 x 10 18 g s" 1 for the whole 
disc. 

However, if the disc is irradiated as mentioned earlier, 
larger disc radius may have a non-negligible contribution to 
the U-band brightness. We therefore use two values (10 1 
and 10 12 cm) for the outer disc radius and study the differ- 
ences. 



5.2 Histogram 

We have analyzed our data in two different ways, firstly us- 
ing cummulative histograms and secondly using power den- 
sity spectra. We use the U-band data so that we can compare 
both methods directly. In this section we compute the cu- 
mulative histogram of the flare duration t from our observed 
data. We define the flare duration as the distance between 
two adjacent minima of a flare. For a flare to be selected it 
must be sampled by a minimum number of points. Because 
we smooth the data using a 3-point box-car function, hence 
this limit must be larger than 3. However, if we set this limit 
to be too large, short flares will be missed. After a detailed 
investigation, we decided to take a minimum of 3 points ei- 
ther side of the flare peak. Therefore, the shortest flare will 
contain 7 points and, because the time resolution is ~45 s, 



Searching for flickering statistics in T CrB 5 




log r [cm] 

Figure 2. Critical mass accretion rates for the accretion disc in 
T CrB as a function of the distance from the central star r. 
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Figure 3. Normalized cumulative histogram of the observed 
flare duration, together with 6 simulated histograms for different 
a values; from bottom to top - 0.001, 0.01, 0.1, 0.5, 1 and 10. 
The simulated histograms have been computed for the following 
parameters: Mi = 1.37 M Q , M acc = 1.9 X 10 18 g s _1 , R out = 
10 11 cm and R in = 3.5 X 10 9 cm. 

they will last for ~315s. In Fig. [3] we present a cumulative 
histogram of the flares, from 300 s onwards. 

In our simulations we need to choose a time step At to 
calculate the final distribution function of duration t. Using 
At = 200 s and a minimum flare duration of t = 300 s we 
obtain a synthetic cumulative histogram which can be com- 
pared with the observed data using a Kolmogorov-Smirnov 
test (KS test). This test searches for the maximal difference 
D between two cumulative histograms; 




-3 -2.5 -2 -1.5 -1 -0.5 0.5 1 

log(alpha) 

Figure 4. Contour plots of KS test for KS values 0.3, 0.2 and 0.1 
(from left to right). 

D = max|Oi(£) - S\{t)\ (23) 

where Oi(i) is the observed histogram and Si(t) is the sim- 
ulated distribution function/histogram. 

Some examples of simulated histograms are shown in 
Fig. |3] The simulations were computed for the parameter a 
in the interval 0.001 - 10. Fig. [4] shows the results of the 
KS test in the a-7?i n plane. It is clear that this method is 
not sensitive enough to constrain the two free parameters. 
The best KS test value 0.1 is too wide and so there are no 
closed contours. Reducing the KS value on the Fig. [3] down 
to ~ 0.0001 does not show more contours features and so 
this method is clearly is not suitable for our study. 

5.3 Power Density Spectrum 

Our method of flare identification can be very subjective 
because different timescales are superimposed. Therefore we 
decided to analyze our data in a another way by calculat- 
ing the power density spectra (PDS) of our observations 
and the synthetic data. We applied the Lomb-Scargle al- 
gorithm (Scargle 1982) which is adequate to deal with un- 
evenly spaced data. Fig. [5]show PDS of our data in UVBR 
Johnson filters. The PDS of the B and V-band lightcurves 
are very similar. Whereas the PDS for the other filters show 
some fluctuations which can be caused by the larger scat- 
ter/noise in U color and the lower amplitude of the vari- 
ability in the J?-band. Therefore, we decided to focus on the 
V-band data, which is also adequate after the outer disc 
estimate in Sec. 15.11 

Fig. [7] shows examples of PDS for different viscosity 
and outer disc radius. The value for the outer disc radius 
of interest are 10 11 and 10 12 cm. Differences in R ou t only 
start to be visible for high value of a. The two R ou t cases 
are indistinguishable for a — 0.5 and 2.5 (used as upper 
limit of viscosity in the following PDS study). Only for a = 
5.0 can the two radii be distinguished. We therefore adopt 
this viscosity value as our sensitivity limit, for the outer 
disc radius in our interval of interest. According to the disc 
instability model (see e.g. Lasota 2001 for review) accretion 
discs typically have a < 1.0, hence our interval of interest 
is lower than the sensitivity limit. Therefore, this parameter 
has no influence in our simulation, and so we fixed it to 
-Rout = 10 11 cm. 



6 A. Dobrotka et al. 



V filter 
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Figure 5. The power density spectrum of our observed data in 
different Johnson filters compared to the V-band data. 

For a given pair of the input parameters 7?i n and a we 
produced 100 synthetic runs with the flare duration statis- 
tics calculated using the method described in Sec. U The 
flares has a triangular shape and the time sampling of the 
runs were sampled the same as the observed runs (a few ex- 
amples are depicted in Fig. 1 10p . We then calculated a mean 
PDS. Fig. [6] shows six examples of the simulated PDS for 
a — 0.001 — 10. We subsequently calculated the difference 
in x 2 between the observed and simulated PDS as follows 

N t \2 

X^E 1 ^ (24) 

i=l 

where Oi and d are the observed and calculated values re- 
spectively, ei is the error of the observed PDS value and N 
is the number of bins in the PDS (we used 15). We then 
computed a 2-D \ 2 grid fits of the inner disc radius and vis- 
cosity parameter (Fig. [8J for the same system parameters as 
for the KS test. This method is clearly more suitable for our 
study as it is more sensitive to the change in input param- 
eters and the wide plateau seen the KS test is not present. 
On the contrary the contour-features are well localized and 
show that the input parameters a and R ln are correlated at 
some interval. 

A minimum \ 2 of 282 is found at Ri n — 3.9 x 10 9 cm and 
a ~ 0.9. Fig.[8]shows the 90% probability (1.6-ct) confidence 
levels in the a — Ri n plane after rescaling so that the reduced 
X 2 is 1. The parameters show a linear trend, an indication 
that the two are correlated. The correlation is constrained 
to the interval i?i n — 1.1 — 4.3 x 10 9 cm for all acceptables 
a values for a disc in hot branch, i.e. 0.1 - 1.0. 



6 DISCUSSION 

The 90% confidence level contours of the PDS fits demon- 
strates that the inner disc radius and the viscosity parameter 
are correlated in the interval -Rm ~ 1.1 — 4.3 x 10 9 cm for 
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Figure 6. The mean power density spectrum of the observed 
runs (the solid points with error bars) compared with 6 ex- 
amples of simulated power density spectra (solid line) for a = 
0.001, 0.01, 0.1, 0.5, 1 and 10 (from bottom to top). 




a= 0.50 




-3.5 -3 -2.5 -2 

log(f) [Hz] 

Figure 7. Examples of simulated power density spectra for dif- 
ferent viscosity and outer disc radii. The inner disc radius has 
been fixed to i?; n = 3.5 x 10 9 cm. 



all aceptable a values for a disc in hot branch, i.e. 0.1 - 1.0. 
Our best solution Ri n — 3.9 x 10 9 cm and a ~ 0.9 is in agree- 
ment with the inner disc radius obtained by Skopal (2005) 
through modeling of UV/optical/IR spectra of TCrB; i.e. 
0.05 - 0.10 R© ~ 3.5 - 7.0 x 10 9 cm. This is larger than 
the estimated white dwarf radius and hence the disc seems 
to be truncated. Furthermore, the best model viscosity pa- 
rameter is in agreement with expectations for a disc in the 
permanent hot branch, i.e. 0.1 < a < 1.0. 
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log(alpha) 

Figure 8. The \ 2 contours of the power density spectra at the 
90% confidence (1.6-c) level. 



Our best values are derived from the minimal value of 
X 2 • The 90% contours levels suggest a large uncertainty in 
parameters which are correlated. The determination of one 
of the two parameters requires an independent measurement 
of the other. Taking the inner disc radius from Skopal (2005) 
to lie in the interval of 3.5 < Ri n < 4.3 x 10 9 cm then gives a 
viscosity interval of 0.7 < a < 1.0. King et al. (2007) report, 
for fully ionized geometrically thin accretion discs, that the 
observational evidence suggests a ~ 0.1 — 0.4. In particu- 
lar, the simulations for SS Cyg (Schreiber et al. 2003) yield 
a viscosity value of about 0.1 — 0.2, whereas Schreiber and 
Gansicke (2002) obtain a = 0.5 for SS Cyg in the ionized hot 
state. Our derived viscosity value of 0.7— 1.0 is considerably 
higher than previous values. However, Schreiber et al. (2003) 
used a maximum mass transfer rate of ~ 1 x 10 17 g s _1 in his 
simulations. This is one order of magnitude lower than mass 
accretion rate in T CrB. The high mass transfer rate through 
the disc together with the irradiation by the central source 
may explain the higher temperature and viscosity than in 
standard dwarf nova system. Furthermore the lower viscos- 
ity values were determined from the long-term light curve 
behaviour where the considered timescales are typical for 
the largest disc radii, whereas our constraints apply mainly 
at the inner disc radius. This gives an idea that the viscos- 
ity parameter a can also be a function of the distance r. In 
the standard disc instability model (see e.g. Lasota 2001 for 
review) the main difference between the viscous parameter 
in the cold and hot branch is due to the ionisation of hy- 
drogen. Once the material is in the hot state, the ionization 
state of hydrogen does not change with increasing tempera- 
ture and the viscosity parameter is supposed to be constant 
in the whole disc. Intensive studies of the disc behaviour 
with constant a were performed in the case of dwarf novae 
and low mass X-ray binaries. But the typical disc radii of 
these binaries are of the order of 10 10 - 10 11 cm. The case of 
symbiotic systems, such as T CrB is different because of the 
large binary separation and the possible outer disc radii up 
to 10 12 cm. But we showed that our simulations are not sen- 
sible to such large disc radii in the acceptable a regime. The 
outer disc radius used for T CrB is 10 11 cm, and hence is in 
agreement with the radius of low mass X-ray binaries mod- 
elled by Lasota (2001) with constant a values. Therefore we 
conclude that radial variations of the viscous parameter is 
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Figure 9. A simulation of the scaling parameter k, disc scale 
height H and radial velocity v x for Mi = 1.37 Mq, M acc = 1.9 X 
10 18 g s" 1 and a = 0.5. 

not required. An alternative explanation for the different a 
results could be a different model other than the Shakura 
and Shunyaev (1973) model for the large outer discs in sym- 
biotic systems. An application of our simulations to dwarf 
novae with disc radii of 10 10 cm can test whether the a differ- 
ences remain or disappear. This case is problematic because 
dwarf novae discs are not in a steady-state. Another pos- 
sibility is to study dwarf novae in outburst when the disc 
is fully ionised and is in a quasi steady-state (see e.g. La- 
sota 2001) and in a permanent hot state except for in the 
lightcurve minima e.g. in VYScl systems (see e.g. Warner 
1995). 

Another interesting question is the location of the flick- 
ering events in the inner disc. This can be understood from 
the behavior of the scaling parameter k(r) (Fig. [9]). For very 
small r values, in the vicinity of the central star, k(r) rises 
very steeply because of the small disc scale height (Fig. [9j 
and hence small turbulent eddies. Furthermore, k(r) de- 
creases with increasing r and the disc scale height, and thus 
less events are needed to transport the angular momentum. 
This strongly suggests that the inner parts of the disc are 
the source of the fast flickering, which is in agreement with 
observations (see e.g. Bruch 1992, Zamanov & Bruch 1998, 
Bruch 1996, Bruch 2000, Baptista & Bortoletto 2004, Za- 
manov 2004). A similar behavior is deduced from the sen- 
sitivity of our modeling to the input parameters R[ n and 
Rout, which is very high for Ri n ~ 10 9 — 10 10 cm but not for 
Rout ~ 10 11 - 10 12 cm. This suggests that the majority of 
turbulent events are located in the central part of the disc 
up to radial distance ~ 10 10 cm. 

The inner disc radius and the viscosity is correlated 
only for low values of r comparing to the disc radius. This 
can be explained because the viscosity must increase with 
higher Ri n in order to get the same \ 2 residual. This means 
that by truncating the central part of the disc, we are re- 
moving many fast events which should be compensated by 
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the increase of the radial velocity of slower events at larger 
radii. The radial velocity v T ~ a 4 ^ 5 strongly decreases at 
very small r (see Fig. [9} . An alternative explanation relies 
on the total number of events in the central disc. The scaling 
parameter k gives an estimate of the number of events and 
this is very large at small r. Hence by truncating the central 
part of the disc we are removing many events. Whereas at 
larger r the decrease of k and v r is flatter and the correlation 
between Ri n and a is not as strong as for small values of r. 

Fig. [TO] shows three examples of simulated lightcurves 
computed for different viscosities. It is obvious that high fre- 
quency fluctuations increase with a, whereas low- frequency 
fluctuations decrease with a. This disappearance of large 
events is visible in our PDS (Fig. [6} ■ The low viscosity PDS 
has a strong peak near frequency 10 -3 Hz which is caused by 
the quasi-periodic oscillations present in the simulated data. 
Such variations has a timescale of ~ 0.01 day clearly visible 
on Fig.JTUJfor a = 0.01 case. The timescale of the data is too 
short to distinguish whether the fluctuations are stochastic 
or quasi-periodic. For higher viscosities, large events disap- 
pear and so does the evidence for quasi-periodic signals in 
the simulated data. Hence the peak at 10 -3 Hz dissapear. 

In addition, the observed PDS in Fig. [6] are scattered 
which produce a high % 2 - This behavior tells us something 
about the observational demands. Our observed PDS is com- 
puted from 9 runs, whereas the simulated PDS is computed 
from about 100 synthetic runs. Smoother observed PDS 
(made from a higher number of observational data sets) are 
required for the analysis to get lower x 2 and then make the 
method more sensitive to the input parameters. The mean 
duration of our data is ~ 2.6 hours. Longer observations 
could remove the quasi-periodic signals at low viscosities 
and reduce the scatter in the observed and simulated PDS 
at low frequencies (~ 10 -4 Hz). Using of a single filter can 
improve the data resolution to ~ 10 seconds and hence the 
quality of the PDS at higher frequencies (~ 10~ 2 Hz). Also 
the observational gaps needed to observe a comparison star 
(clearly visible on Figs. [T1and ll0[) affect the PDS. Therefore 
CCD observations with fast readout speeds are preferred. 

Our best model yields a parameters in reasonable agree- 
ment with other works. We therefore conclude that our sim- 
ple idea of angular momentum transport through turbulent 
eddies is sufficient to explain the observed flickering statis- 
tics. No other mechanisms are then required to generate the 
observed nickering activity. 



7 SUMMARY 

We have analyzed U-band photometry of the aperiodic vari- 
ability of T CrB during the high state, taken using data with 
a time resolution of ~ 45 s. We developed a method to sim- 
ulate flickering activity based on a simple idea of angular 
momentum transport through turbulent processes in the ac- 
cretion disc. We simulated artificial data sets and compared 
them with the observations using cumulative histogram with 
Kolmogorov-Smirnov tests and power density spectra. The 
best model yielded Tim — 4 x 10 9 cm and a ~ 0.9 but with 
large uncertainties. However, using the results from litera- 
ture for the inner disck radius, we derived the viscosity to 
lie in the interval 0.7 < a < 1.0. The disc is weaklu trun- 
cated and we provide an explanation for the concentration 
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Figure 10. Three examples of the synthetic data compared to 
an observed lightcurve. The simulated runs are computed for Mi 
= 1.37 M Q , Ma.cc = 1.9 X 10 18 g s" 1 and R in = 3.5 X 10 9 cm and 
three different viscosities a = 0.01, 0.50 and 10.0. 

of flickering events in the central part of the accretion disc. 
The correlation between the parameters Ri n and a is ex- 
plained by the coupling between the number of fast events 
and the radial velocity, i.e. the removal of fast events pro- 
duced by increasing R[ n is compensated by an increase in 
v T and hence a. Our turbulent scenario of angular momen- 
tum transport successfully explains the observed flickering 
statistics without any other mechanism required. 
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